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ABSTRACT 

We present physical properties of two submillimeter selected gravitationally lensed sources, identi- 
fied in the Herschel Astrophysical Terahertz Large Area Survey. These submillimeter galaxies (SMGs) 
have flux densities > 100 mJy at 500 /im, but are not visible in existing optical imaging. We fit light 
profiles to each component of the lensing systems in Spitzer IRAC 3.6 and 4.5 jam data and success- 
fully disentangle the foreground lens from the background source in each case, providing important 
constraints on the spectral energy distributions (SEDs) of the background SMG at rest-frame optical- 
near-infrared wavelengths. The SED fits show that these two SMGs have high dust obscuration with 
Ay ~ 4— 5 and star formation rates of ~ 100 M yr _1 . They have low gas fractions and low dynamical 
masses compared to 850 /im selected galaxies. 

Subject headings: galaxies: individual (SDP.81: H- ATLAS J09031 1.6+003906, SDP.130: H- ATLAS 
J091305. 0-005343) — galaxies: starburst — gravitational lensing: strong 
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1. INTRODUCTION 

Gravitational lensing is an invaluable astrophysical 
tool. It can be exploited to study galaxies beyond instru- 
mental blank field sensitivities and to constrain the total 
mass of galaxy systems, without regard for its dark or 
luminous nature (see review by Treu 2010). The magni- 
fication of distant galaxies through gravitational lensing 
enables the detailed study of sources that would other- 
wise be undetectable. In the submillimeter regime this 
includes members of the population of intrinsically faint 
galaxies that comprise a significant fraction of the cosmic 
far-infrared background. 

The Herschel 35 (Pilbratt et al. 2010) Astrophysical 
Large Area Survey (H-ATLAS; Eales et al. 2010) is the 
largest open-time key project currently being undertaken 
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by Herschel and aims to survey 550 deg 2 and detect 
> 300,000 galaxies. A new methodology for selecting 
gravitational lenses using wide-area submillimeter (sub- 
mm) surveys (Blain 1996, Perrotta et al. 2002; Negrello 
et al. 2007) has been tested for the first time during Her- 
schel^ science demonstration phase (SDP; Negrello et 
al. 2010, N10 hereafter). Candidates are selected us- 
ing a flux cut of 100 mJy at 500 /xm, a limit based on 
the steep number counts slope in the sub-mm (Negrello 
et al. 2007). Above this flux limit only gravitationally 
lensed objects and easily identifiable 'contaminants' re- 
main e.g. blazars and local spiral galaxies. Initial results 
show that this selection method has ~100% efficiency 
and should deliver a sample of hundreds of new grav- 
itational lenses in planned wide area Herschel surveys, 
probing galaxies with intrinsic fluxes below the Herschel 
confusion limit. 

Five lens candidates in the SDP H- ATLAS data (Pas- 
cale et al. 2010; Rigby et al. 2010; N10) were confirmed 
with spectroscopic redshifts obtained via the detection of 
carbon monoxide (CO) emission lines (Lupu et al. 2010; 
Frayer et al. 2011) of the background galaxies and opti- 
cal spectra (Negrello et al. 2010) of the lens galaxies. In 
this Letter we study two of these gravitational lenses, 
H-ATLAS J090311.6+003906 (SDP.81) and H-ATLAS 
J091305.0-005343 (SDP.130), which are submillimeter 
galaxy (SMG) at redshifts of z = 3.04 and z = 2.63 
being lensed by ellipticals at z — 0.30 and z — 0.22, 
respectively. Submillimeter Array (SMA) imaging re- 
veals the sub-mm morphology, consistent with a lensing 
event, with multiple peaks distributed around the po- 
sition of the foreground elliptical galaxy (N10). While 
no lensed background images were detectable in optical 
imaging, the spectral energy distribution (SED) models 
suggest that flux from the background sources should be 
detectable above ~10 /iJy at near-IR wavelengths, from 2 
to 5 fim (N10). At these wavelengths, the emission from 
the foreground lenses, at z ~ 0.25, and the background 
SMGs becomes comparable. 

In this Letter we present Spitzer /IR AC imaging, light 
profile models, photometry and physical characteristics 
for SDP.81 and SDP.130 derived from SED fitting. In the 
next section, we present a summary of the Spitzer data, 
while in Section 3 we discuss modeling of the light pro- 
file of the foreground lenses. We perform SED modeling 
of background SMGs and present results related to the 
properties of these two sources in Section 4. Through- 
out the Letter we assume flat ACDM cosmology with 
^ M = 0.3 and H = 70 kms" 1 Mpc -1 . 

2. SPITZER IR AND KECK OPTICAL IMAGING DATA 

The IR imaging data are part of Spitzer program 548 
(PI: A. Cooray), released for analysis on 2010 July 30. 
For this program, Infrared Array Camera (IRAC; Fazio 
et al. 2004) images of SDP.81 and SDP.130 were taken 
at 3.6 /im (Channel 1) and 4.5 /am (Channel 2). Both 
lensing systems were imaged with a 36-position dither 
pattern and an exposure time of 30 s for each frame to 
achieve an effective total exposure time of 1080 s. The 
rms depth reached is 0.3 and 0.4 /aJy in Channels 1 and 

35 }{ ersc }i e i i s an ESA space observatory with science instru- 
ments provided by European-led Principal Investigator consortia 
and with important participation from NASA. 
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Fig. 1. — Keck i-band (top row), IRAC Channel 1 (middle row), 
and Channel 2 (bottom row) postage stamp images for SDP.81. 
The corresponding residual (right column) is presented after sub- 
traction of single de Vaucouleurs profiles modeled on the Keck i- 
band data. The Keck residuals show no signs of structure asso- 
ciable with the overlaid 880 fim SMA signal-to- noise ratio (S/N) 
contours. For the IRAC residual there is structure that can be 
associated with the SMA contours, which suggests this is lensed 
structure of the respective background galaxy. The SMA contours 
are overlaid in steps of 3cr, 8cr, 13a etc. SMA peaks with S/N > 
8 (four for ID81 and two for ID130) have a maximum rms on the 
position of a point source of 0.1". All postage stamps are plotted 
with the same pixel intensity scaling after normalizing the Keck 
z-band data to the dynamic pixel value range of the IRAC data. 



Corrected basic calibrated data pre-processed by the 
Spitzer Science Center (SSC), using the standard pipeline 
version S18.18.0, were spatially aligned, resampled, and 
combined into a mosaic image using version 18.3.1 of 
the SSC's MOPEX software suite (Makovoz & Marleau 
2005). The IRAC mosaics have a resampled pixel size of 
0.6" and angular resolution of 2 to 2.5". For the work 
presented here, we also use the IRAC point spread func- 
tion (PSF; version 2010 April) file as provided by the 
SSC 36 . 

We also make use of optical images of SDP.81 and 
SDP.130 that were acquired on 2010 March 10 using the 
dual- arm Low Resolution Imaging Spectrometer (LRIS; 
Oke et al. 1995, McCarthy et al. 1998) on the Keck I 
telescope and reported in N10 (see Figures [1] and [2]). 
Each target received simultaneous 3 x 110 s integrations 
with the ^-filter and 3 x 60 integrations with the z-filter 
using the blue and red arms of LRIS, respectively. A 

36 http : //ssc . spitzer . caltech . edu/irac/calibrationf iles/psf prf / 
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20" dither pattern was employed to generate on-sky flat- 
field frames. We performed photometric calibration us- 
ing 1 s g- and z-band observations of bright stars in each 
field. The data were reduced using IDL routines and 
combined and analyzed using standard IRAF tasks; the 
seeing FWHM of the final science exposures is ~0.8". 
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Fig. 2. — Keck z-band (top row), IRAC Channel 1 (middle row) 
and Channel 2 (bottom row) postage stamp images for SDP.130. 
The corresponding residuals are presented after subtraction of sin- 
gle de Vaucouleurs profiles modeled on the Keck i-band data. Al- 
though the Keck residuals show no signs of significant structure, 
the subtracted IRAC data have residual structure that can be as- 
sociated with the SMA contours, which suggests lensed structure 
of the respective background galaxy is present. The SMA contours 
are presented as in Figure [T] All postage stamps are plotted with 
the same pixel intensity scaling after normalizing the Keck z-band 
data to the dynamic pixel value range of the IRAC data. 



3. MODELING THE LENSES 

To construct models of the light profiles for each lens- 
ing system we use GALFIT (Peng et al. 2002), which 
allows multiple profiles per object and performs a si- 
multaneous nonlinear minimization. Prior to fitting pro- 
files to SDP.81 and SDP.130 in the IRAC data, de Vau- 
couleurs profile (de Vaucouleurs 1948) models are con- 
structed for the Keck i-band imaging, to take advantage 
of the comparatively higher resolution of these images, 
and the presence of only one component per system i.e. 
the lens galaxy. 

To look for any potentially lensed structure in the 
IRAC data, the Keck models are used to fit the IRAC 
Channels 1 and 2 data, keeping the effective radius and 
ellipticity fixed, and using the appropriate IRAC PSF 
36 . On subtraction of the results, and in comparison to 
the model subtracted Keck data and SMA contours, the 
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Fig. 3. — Multi-component light profile models for SDP.81 (top) 
and SDP.130 (bottom), for the IRAC Channel 1 data. Shown are 
the un-subtracted data, the lens model (a single Sersic profile for 
SDP.81 and a Sersic plus exponential disk profile for SDP.130 ), 
the residual after subtraction of the lens model profile only and 
the total model residual (full residual). The SMA contours are 
presented as in Figure [T] All data are plotted with the same pixel 
intensity scaling. 

IRAC band residuals strongly suggest that a more com- 
plex structure, associated with the background SMG, is 
present for both SDP.81 and SDP.130 (see Figures [1] and 

ED. 

We verify that these residuals are significant, and not 
an artifact of imperfections in the IRAC PSF, by compar- 
ing them with residuals derived for three (non-lensing) 
elliptical-like galaxies in the same field, after fitting them 
with single Sersic profiles. For both the lensing systems 
and the comparison ellipticals aperture flux ratios were 
determined for the residual image and the corresponding 
un-subtracted data. To consider only positive structure, 
pixel values > 2 a below the local background were re- 
placed with the median local sky value. The SMG resid- 
uals were found to have flux ratios ~ 3 — 5 greater than 
those for the random elliptical galaxies. 

In order to more precisely disentangle the lens and 
background components we represent the flux from the 
lensed SMGs with Sersic profiles. Peaks were identified 
in the IRAC single profile fit residual images, by fitting 
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Gaussian profiles. We add three profiles for SDP.81 and 
two for SDP130, which all correspond to significant sub- 
mm contour peaks in the SMA data, within 0.8". N10 
modeled the Keck % band with a Sersic profile plus an 
exponential disk, as this combination gave a marginally 
improved \ 2 over single profile models, and found that 
the exponential component significantly contributes to 
the total profile for SDP.130. Therefore, a Sersic plus ex- 
ponential disk profile is adopted for the SDP.130 lens and 
a single Sersic profile for the SDP.81 lens. To produce 
the final models for the lens and background galaxies, 
we refit the profile representing the lens galaxy and the 
additional Sersic profiles simultaneously for each system 
with GALFIT. These final model fits for the background 
component of SDP.81 show a partial "Einstein ring" -like 
morphological structure for the lensed source (Figure. EJ). 
For SDP.130, the lensed component is more compact and 
in close proximity to the lens profile. The resulting best 
fit profiles for the background galaxies agree well with 
the SMA contours, and the combined models subtract 
cleanly, suggesting successful lens/source decoupling. 

4. SPECTRAL ENERGY DISTRIBUTION OF THE LENSED 

SMGS 

The GALFIT-integrated magnitudes for each com- 
ponent of the final SDP.81 and SDP.130 models were 
converted to flux density to extend the existing multi- 
waveband photometry (see table [T] and N10, and refer- 
ences therein) into the near-IR. We assign 1 a errors to 
the photometry obtained from the final GALFIT model 
profiles, using the magnitude distributions for all the 
GALFIT trials that converged. PACS re-imaging of the 
lensed H-ATLAS sources provides new photometry at 
160 /im and upper limits at 70 fim (I. Valtchanov et al. in 
preparation; Ibar et al. 2010). For the goal of deriving 
physical properties, the IRAC photometry adds partic- 
ularly important constraints to the SMG SEDs, which 
previously consisted of just upper limits at wavelengths 
below 250 /mi for SDP.130 and 160 /im for SDP.81. 

The SEDs of the SMGs are fitted using the models 
of da Cunha et al. (2008), calibrated to reproduce the 
ultraviolet-to-infrared SEDs of local, purely star-forming 
Ultra Luminous Infrared Galaxies (ULIRGs; 10 12 < 
L m /L e < !0 13 ; da Cunha et al. 2010). The SED 
models assume a Chabrier (2003) initial mass function 
(IMF) that is cutoff below 0.1 and above 100 M ; using 
a Salpeter IMF instead gives stellar masses that are a 
factor of ~ 1.5 larger. We find that a significant atten- 
uation by dust (Ay ~4-5) is required to be consistent 
with the IRAC photometry and optical/near-IR upper 
limits (Figure EJ, which is consistent with other ULIRGs 
and SMGs (e.g. Geach et al. 2007; Hainline et al. 2010; 
Michalowski et al. 2010; Wardlow et al. 2010). 

Using the Chabrier (2003) IMF, with parameters de- 
rived from the SED fits, we find that SDP.81 and 
SDP.130 have stellar masses (M*) of (2.5±1.7) x 10 11 M© 
and (4.5±2.5)xlO n M , respectively (see tabled}. How- 
ever, we note that there is an additional systematic error 
of up to a factor of 10, due to uncertainty in the appro- 
priate mass-to-light ratio (see Wardlow et al. 2010 for a 
discussion) and magnification factors, as well as a lack of 
observations at optical/near-IR wavelengths. J and Ks 
photometric data on these galaxies with VLT/HAWK-I 
observations (A. Verma et al. in preparation 2011) could 



TABLE 1 

Parameters of the Foreground Lens and Background SMG 
for SDP.81 and SDP.130 



Parameters 


SDP.81 


SDP.130 


Foreground lens 


RA 


09 h 03 m 11.6 s 


09 h 13 m 05.0 s 


Dec 


00 d 39 m 06 s 


-00 d 53 m 43 s 


a 

Redshift 


0.299 


0.220 


SDSS u (/xJy) 


3.9 ± 2.0 


1.7 ± 1.7 


SDSS g (/iJy) 


24.9 ± 1.1 


19.4 ± 0.7 


SDSS r (/xJy) 


115 ± 2 


66.1 ± 1.2 


SDSS i (/xJy) 


198 ± 4 


109 ± 2 


SDSS z (/xJy) 


278 ± 8 


143 ± 7 


UKIDSS Y (/xJy) 


320 ± 20 




UKIDSS J (>Jy) 


370 ± 20 




UKIDSS H (>Jy) 


510 ± 50 




UKIDSS K (/iJy) 


570 ± 70 




Spitzer 3.6 /im (mJy) 


0.35 ± 0.04 


0.213 ± 0.03 


Spitzer 4.5 /im (mJy) 


0.22 ± 0.04 


0.230 ± 0.01 


Background SMG: observed quantities 


a, 

Redshift 


3.042 


2.625 


a 

Magnification 


25 ± 7 


6 ± 1 


Keck/LRIS g (/iJy) 


< 0.13 


< 0.20 


Keck/LRIS i (/iJy) 


< 0.20 


< 0.35 


UKIDSS Y (/xJy) 


< 6.27 




UKIDSS J (/iJy) 


< 9.23 




UKIDSS H (/iJy) 


< 8.52 




UKIDSS K (/iJy) 


< 13.5 




Spitzer 3.6 /im (mJy) 


0.062 ± 0.04 


0.044 ± 0.01 


Spitzer 4.5 /im (mJy) 


0.126 ± 0.05 


0.047 ± 0.01 


PACS 70 /im (mJy) b 


< 8.0 


< 9.0 


PACS 100 /im (mJy) 


< 62 




PACS 160 /im (mJy) b 


51 ± 5 


45 ± 8 


SPIRE 250 /im (mJy) 


130 ± 20 


105 ± 17 


SPIRE 350 /im (mJy) 


180 ± 30 


128 ± 20 


SPIRE 500 /im (mJy) 


170 ± 30 


108 ± 18 


SMA 880 /im (mJy) 


76 ± 4 


39 ± 2 


IRAM 1200 /im (mJy) 


19.6 ± 0.9 


11.2 ± 1.2 


Background SMG: derived quantities 


L m (lO 12 L )^ 


2.0 ± 0.6 


5.6 ± 1.2 


A v 


4.4 ± 0.6 


5.0 ± 0.5 


M H 


-25.7 ± 0.3 


-26.01 ± 0.17 


SFR (M yr- 1 ) 


74 ± 30 


150 ± 50 


M(H 2 ) (10 10 M ) d 


1.4 


2.7 


M* (10 11 M ) e 


2.5 ± 1.7 


4.5 ± 2.5 


M dust (10 8 M ) f 


3.4 ± 1.0 


11 ± 2 


g 

fi 


0.05 ± 0.01 


0.08 ± 0.01 



Note. — The optical and sub-mm flux densities and 3a upper lim- 
its are taken from N10, unless otherwise noted. Derived quantities are 
corrected for magnification and the errors quoted account for the tab- 
ulated uncertainty in the magnification. 

(a) Spectroscopic redshifts (N10). 

(b) From PACS re-imaging data (I. Valtchanov et al. in preparation 
2011; Ibar et al. 2010). 

(c) 2 to 1000 /im luminosity based on the best-fit SED; 

(d) Frayer et al. 2011, uncertainty is at least a factor of ~ 2. 

(e) There is an additional systematic uncertainty up to a factor of ~ 10 
(see text). 

(f) Calculated with a dust mass absorption coefficient approximated as 
K\ oc with (3 — 1.5 for warm dust and (3 = 2 for cold dust, and 
a normalization of 0.77 g _1 cm 2 at 850 /im. There is an additional 
systematic uncertainty of at least a factor of ~ 2. 

(g) Gas fraction, /i = M(Hq,) /[M* + M(H2)], uncertainty is a factor of 
2-3. 
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Fig. 4. — Left: photometry and best-fit SEDs for the foreground elliptical (blue) and background SMG (red) for SDP.81 and SDP.130. 
The photometric points and upper limits are taken from Negrello et al. (2010), with updated PACS flux densities at 160 /im and upper 
limits at 70 /mi. Photometry from the best fitting light profiles to the IRAC data are added for the foreground lens galaxies (turquoise 
points) and for the background SMG (pink points). The SEDs are fitted using the models of da Cunha et al. (2008). We find that high levels 
of visual extinction of the SMGs are required (Ay > 4) to be consistent with the optical-to-NIR data. Top right: a plot of (magnification 
corrected) rest-frame absolute if-band (Mh) magnitude against redshift for SDP.81 and SDP.130, showing the approximate correspondence 
of Mh with stellar mass. For comparison, we also show 850- and 870-/im selected galaxies (Hainline et al. 2010; Wardlow et al. 2010). The 
two lensed H- ATLAS sources are brighter than ~ 95% of the comparison sample, and correspondingly, are likely to be amongst the most 
massive. Error bars are omitted from the comparison sample for clarity. Bottom right: ratio of CO and if -band luminosity against CO 
line velocity for the H- ATLAS lensed galaxies and 850 /mi selected SMGs. Lco/^H is independent of the lensing model (to the first order; 
see discussion in the text) and represents the gas-to-stellar mass ratio, crco is indicative of dynamical mass, although it also depends on 
the inclination angle and size of the CO emitting region, which increases the observed scatter. SDP.81 and SDP.130 have small dynamical 
masses and small gas fractions, relative to the comparison sample. 



potentially improve the estimates of extinction and stel- 
lar mass. 

Rest-frame absolute iLband magnitudes (Mh) provide 
a guide to galaxy stellar masses that is not complicated 
by the details of the assumed star formation history and 
is straightforward to compare to other similarly selected 
samples. In the top panel of Figure 3] we plot Mh against 
redshift for SDP.81 and SDP.130 compared to 850 and 
870 /im selected SMGs (Hainline et al. 2010; Wardlow 
et al. 2010). The H- ATLAS lensed galaxies are brighter 
than the majority of the 850-870 /im selected SMGs, 
thus if there is no active galactic nucleus contribution to 
their iL-band luminosities they are likely to be amongst 
the most massive. 

We next consider the CO(l— 0)-to-ff-band luminosity 
ratio (Leo/ Lb) , a representation of the gas-to-stellar 
mass fraction that is mostly independent of the lensing 
model. However, the IRAC data trace stellar emission, 
which may originate from a spatially different region of 
the SMG to the far-infrared emission. As such there 
may be variation in the lensing amplifications across the 
source plane, potentially leading to a small effect on 
Leo /Ln The bottom panel of Figure |4] shows Lqo/Lyl 
against the CO linewidth (crco) f° r the two H- ATLAS 
lensed galaxies, compared to 850-/im selected SMGs 



(Hainline et al. 2010). crco is independent of Lqo/Lh 
and is indicative of the dynamical mass of the system, al- 
though some scatter is introduced due to the dependence 
on the inclination angle and the size of the CO emitting 
region. Published CO luminosities of the 850 /im selected 
sample (Greve et al. 2005, Coppin et al. 2008; Frayer 
et al. 2008, Tacconi et al. 2008; Bothwell et al. 2010) 
are converted to the equivalent CO (1—0) values using 



J CO(3-2)/LcO(l-0) 



0.6 and L CO (4-3)/£co(i- 



•o) 



0.6 



(Ivison et al. 2010). Galaxies with poorly defined crco 
are excluded. The H-ATLAS lensed galaxies have low 
Lqo/Lh (representative of gas fraction) and low aco 
(representative of dynamical mass) compared to 850- 
/im selected galaxies (Figure 2]). There may also be a 
weak inverse trend between Leo / and <tco for all the 
SMGs, although a larger sample is required for confirma- 
tion. 

Summary: We have studied the background SMGs of 
two sub-mm bright gravitational lenses, identified in the 
H-ATLAS SDP data. The intrinsic sub-mm flux densities 
of these SMG are below the Herschel confusion noise and, 
therefore, undetectable without the fortuitous lensing by 
foreground ellipticals. The full 550 deg 2 survey area of 
H-ATLAS will recover a sample of > 200 lensed SMGs, 
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and Spitzer follow-up observations will enable us to study 
the physical properties of SMGs over a wide range of 
redshift and far-IR luminosities. Studies such as this are 
necessary to further understand the nature of sources 
that contribute to the bulk of the cosmic far-infrared 
background. 
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